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Abstract 

If the geoid and the satellite position are known accurately, satellite altimetry can be used to 
determine the geostrophic velocity of the surface ocean currents. The purpose ot this 
investigation is to simultaneously estimate the sea surface topography, f, the model tor the 
gravity field, and the satellite orbit. Satellite tracking data from fourteen satellites were 
along with Seasat and Geosat altimeter data as well as surface gravity data for the solution 1 he 
estimated model of f compares well at long wavelengths with the hydrographic model or 
Covariance studies show that the geoid is separable from £ up to degree 9, at which point geoid 
error becomes comparable to the signal of C 

1. Introduction 

The determination of the general ocean circulation is one of the more important applications 
future satellite altimeter missions will provide. If the height of the ocean surface relative to the 
geoid can be measured, the absolute geostrophic surface current velocity at a given location can 
be inferred. The Seasat and Geosat altimeter missions measured the height ot the ocean to a 
precision suitable for this purpose [Tapley et a/., 1982; McConathy and Kilgus, 1987]. Errors 
in the gravity field have limited previous determinations of f for two reasons: 1. C is measured 
with respect to the geoid, an equipotential surface defined by the gravity field model, and 2. 
gravity field errors are the primary limitation in the accuracy of the satellite orbit computation. 

The radar altimeter carried by Earth orbiting satellites measures the range between the radar 
antenna and the instantaneous ocean surface. If the geocentric position of the satellite is known, 
the altimeter measurement can be differenced from the computed satellite height to yield the 
ocean surface height at the sub-satellite point. If there were no forces acting except the Earth s 
gravity and centrifugal force, the ocean surface would coincide with an equipotential surface 
referred to as the geoid, N. However, a number of effects cause the height of tne ocean surface 
to deviate from the geoid height. Among these, ocean currents cause deviations with maximum 
amplitudes of about one meter. If an accurate estimate of N is available, the measured height ot 
the ocean surface (corrected for tides and other effects [Tapley et al., 1982]) and the geoid 
height may be differenced to provide an estimate of f. 

Previous maps of ( with wavelengths comparable to a degree and order 6 spherical 
harmonic model (6000 km) have been determined using Seasat altimeter data [Tai and Wunsch, 
1984; Eneelis, 1986], Each of these approaches adopted a mean sea surface computed using 
Seasat altimeter data which was then differenced with a low degree and order geoid jo give an 
estimate of the long wavelength components of f. The error in the gravity field model has been 
a limiting factor in the accuracy achieved in these investigations. 

The results described here were derived as an adjunct to a concentrated effort to determine 
an improved gravity field model for the Topex/Poseidon altimeter mission [Tapley et al , 1987J. 
The use of altimeter data in a gravity field solution, without correcting for f, can allow the 
oceanographic signal due to f to be aliased into the gravity field. The potential tor 
simultaneously improving the geoid and C using satellite altimetry, surface gravity observations, 
and hydrography has been discussed by Wunsch and Gaposchkin [1980]. Since the influence 
of the gravity field on the geoid and the orbit, and the effect of £, are both present in the 
altimeter measurement, a more rigorous treatment of theproblem would involve estimating t e 
parameters which define both models simultaneously. This investigation demonstrates a joint 
solution for £, the gravity field, and the satellite orbit using satellite altimeter and tracking data. 

2. Methods 

Seasat and Geosat tracking data and altimeter data along with tracking data from twelve 
additional satellites and surface gravity data were processed in a simultaneous solution using a 
modified least squares estimator [Tapley et al., 1987]. The estimated p^metere inc ud^ a set 
of spherical harmonic geopotential coefficients, complete to degree ^nd order 50 with selectwl 
higher order terms, the spherical harmonic coefficients of £, complete to 
and ocean tide parameters. Solar radiation pressure, atmospheric drag, and aoPP ler ^cR in g 
Son coordinates were adjusted as individual satellite parameters. Altimeter biases and scale 
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factors for the significant wave height correction were estimated for Seasat and Geosat also. In 
addition to Seasatand Geosat altimeter data, the different types of satellite tracking data included 

ia«pr ranee data doDDler data, optical data, and altimeter crossover data. 

The stationary ocean surface is composed of the sum of N and £• The Earth s gravitational 
potential can be expressed in spherical harmonic form as: 


(/ = 


1+5 £ ( r ffPi m (sin<piCi m cosrnl + Sj^sinmA] 


L £=2 


(i) 


where C Im , S, m are the normalized spherical harmonic coefficients of degree / and order m pi is 
the product of the gravitational constant and the total mass of the earth, r c is the mean equatorial 
SB the earth r tfU are the radial distance, geocentric latitude and longitude Pleasured from 
thezero meridian,^ ( sind» is the normalize! associated Legendre function (m*0) or he 
ncmnalfz^Le gen dre poly nomi al of degree / (m=0), and is. the truncation degree of the 
geopotential. The position vector (x,y,z) of the sub-satellite point, which defines the geoid 
height, is computed oy iteratively solving [Shum, 1983]: 

W{x,y,z) = Wo and VW / 1 VW I =% W 

where W(x y ,z) = U + T, the total potential at the sub-satellite point, T is the rotational 
potential VTn is value of the total potential corresponding to mean sea level, and Uh is ; a unit 
vector normal to the geoid and passing through the satellite. Thus the set of constant 
parameters Ci and Si , i n Equation (1) are common to both the satellite dynamics and the 
geoid definition (2). A Spherical harmonic model was employed to represent the height of C as - 


= ^ [c'lm cosmA + S/ m sin/nA] Pim (sin0) 


(3) 


/ = lm = 0 


i = i rri — v 

where Ci St are the surface spherical harmonic coefficients of the model for £ In this 
constrained to be zero over land areas and in polar regions fl * |> 70 ). Wht e 
the allmieter data with the land constraints do not form an equally spaced data set. the estimate 
for the coefficients retains the general characteristics of a fully orthogonal parameter set. 

Portions of the Seasat and Geosat altimeter data sets were selected to correspond to the time 
periods when the best ground Packing data was available. Using the Preliminary ^ Topex 
Gravity Field (PTGF2) [Tapley et al., 1987], reference orbits c p^E)* t Sd 
da?a and altimeter crossover daia. In all, 30 days 

1978) and 54 days of Geosat altimeter data (1 1/17/86-12/4/86^ 12/21/8o-l/^4/js /) we e 
processed The altimeter data were edited including the removal of data if the water depth was 
Fess than 2000 m (to remove data with large tidal uncertainties) or if the data was over an ocean 
pench or seamount. The data were corrected for the unmodeled short wavelength features of 
the geoid (> degree 50) and compressed to form altimeter measurements 10 seconds apart. 

There^ was no a priori variance on the coefficients off- The a priori variance of the 
geopotential coefficients, v/ 2 , were computed using Kaula's rule [Kaula, 1966] which can 
stated as vf=a 10' 10 l A , where a is a scale factor equal to 0.04. 

3. Results 

To evaluate the satellite altimeter solution for C, a spherical harmonic m^el of f was 
constructed from hydrographic data. The Levitus ocean surface referenced to 2250 dbars is 
sHnnfprl for this purpose f Levitus 1982]. Grid points on land and in the polar regions (| <j> | 
70 9 ) P were set to Fero and a degree and irder 15 spherical harmonic model was fit to the data 
using^ a Minimum Variance estimator. Using the estimated values for the spherical harmonic 
coefficients C can be computed as a function of latitude and longitude using Equation (3). 
Fieure 1 depicts maps of C computed from both the altimeter and Levitus derived spherical 
harmonic coefficiems complete to degree and order 6. The major gyres of the ocean circulation 
JSdJSSSS £h solutions. Up to degree and order 6, the altimeter and hydrographic models 
of C compare well The agreement between the hydrographic and altimetnc results is better than 
in prev^u^altimeter soluSis for £, suggesting a significant 

simultaneously solving for the gravity field, the satellite orbit, and f. Further details ot 

^^The lm^wemeratn^he geoid by computing new mfemnce orbits 

(using the estimated models) and analyzing the residuals from the compressed altimeter data. 
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was found that the altimeter residuals improved from 225 cm (for a gravity field without 
altimeter data) to 30 cm, indicating the geoid model has been substantially improved. 



Figure la. The estimated model of f to degree and order 6 



Figure lb. The Levitus model of £ to degree and order 6 
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4. Covariance analysis 


One benefit of simultaneously estimating both the geopotential coefficients and the 
coefficients of £ is that the resulting covariance matrix will describe not only the individual 
errors in N and £, but also the correlation of these errors. Considerable effort was expended to 
determine the relative weights of the data and the scaling of the covariance matrix. 

By examining the degree variance of the £ and the degree error variance of N ( I able 1) it is 
seen that beyond degree 9, geoid error is larger than the signal of £. The covariance matnx or 
the coefficients of £ and the geopotential can also be used to compute the standard deviation of 
the error in N and £. The error in N is largest over the land areas (where there is no altimeter 
data) and in equatorial regions. For terms up to degree 6, the geoid height is accurate to about 8 
cm (Table 1). The error in £ mimics the geoid error, the largest error occurs near the equator 
and increases with degree. The magnitude of the error in £ is about 10 cm or e P 
degree 6 (Table 1). In addition to their individual errors, the correlation coefficient fortheerror 
in A! and £ can be analyzed to determine if the two quantities are separable. Figure 2 mows ; tl n e 
geographical variation of the correlation coefficient for terms in the wvmanw 
degree 6 in both N and £. This figure verifies that the solution for N and £ are separable, at this 
degree, since the correlation coefficients are less than 0.6 globally. 


144 




Contour Interval =■= 0.05 


Longitude 

JUT U1LCI V tXL — # . 

Figure 2. The correlation coefficient of the errors m £ and N 
5. Conclusions 

Tt has been shown that a global estimate of £ can be obtained from a joint solution for £and 
the eravitv field using multiple satellite tracking data and altimeter data. The simultaneous 
solution for the £ the gravity field, the satellite orbit, and other geophysical parameters provides 
^substantially improve*! estimate for the gravity field and the long-wavelength signal of the 
eenerafocean emulation. This estimate of £ compares favorably at long wavelengths with 
lyd^Sc determinations of {. The accuracy of .Se solution decreases as shorter wavelength 
features are examined, primarily due to the influence of geoid error. 
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